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ABSTRACT 

We present Sunyaev-Zel'dovich Effect (SZE) imaging observations of three dis- 
tant [z > 0.8) and highly X-ray luminous clusters of galaxies, CI J1226. 9+3332, 
CI J0152.7-1357 and MS 1054.4-0321. Two of the clusters, CI J1226.9+3332 and 
CI J0152.7— 1357, were recently discovered in deep ROSAT x-ray images. Their 
high X-ray luminosity suggests that they are massive systems which, if confirmed, 
CN ■ would provide strong constraints on the cosmological parameters of structure for- 

■ mation models. Our Sunyaev-Zel'dovich Effect data provide confirmation that 

they are massive clusters similar to the well studied cluster MS 1054.4— 0321. 
{Sj . Assuming the clusters have the same gas mass fraction as that derived from 

SZE measurements of eighteen known massive clusters, we are able to infer their 
mass and electron temperature from the SZE data. The derived electron tem- 
peratures are 9.8+ 4 ;g, 8.7+ 4 g, and 10.4^|q keV, respectively, and we infer total 
masses of ~ 2 x 10 14 /ij~Q M Q within a radius of 65" (340 kpc) for all three 
clusters. For CI J0152.7— 1357 and MS 1054.4—0321 we find good agreement be- 
tween our SZE derived temperatures and those inferred from X-ray spectroscopy. 
No X-ray derived temperatures are available for CI J1226. 9+3332, and thus the 
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SZE data provide the first confirmation that it is indeed a massive system. The 
demonstrated ability to determine cluster temperatures and masses from SZE 
observations without access to X-ray data illustrates the power of using deep 
SZE surveys to probe the distant universe. 

Subject headings: cosmology: observations — galaxies: clusters: individual (CI J1226. 9+3332, 
CI J0152. 7-1357, MS 1054.4-0321) — Sunyaev-Zel'dovich Effect — cosmic mi- 
crowave background — techniques: interferometric 
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Introduction 



The existence of galaxy clusters at high redshift can place powerful constraints on the 
physical and cosmological parameters of structure formation models (Bahcall & Cen (1992); 
Luppino & Gioia (1995); Oukbir & Blanchard (1997); Donahue et al. (1998); Eke et al. (1998); 
Haiman et al. (2000)). The greatest leverage is provided by the most massive and distant 
clusters (e.g., Viana & Liddle (1996)). Two distant and highly x-ray luminous clusters were 
recently discovered in deep ROSAT x-ray images: C1J1226. 9+3332, a cluster at redshift 
z = 0.89, was discovered in the WARPS survey (Ebeling et al. 2001; Scharf et al. 1997). 
C1J0152.7-1357, at redshift z = 0.83, was detected in the RDCS, SHARC, and WARPS 
surveys (Delia Ceca et al. 2000; Romer et al. 2000; Ebeling et al. 2000). Based on their 
x-ray luminosities (Lx[0.5 — 2keV] > 2 x 10 44 hJ~ 2 erg s _1 ), these clusters are thought to 
be highly massive which, if confirmed, will provide significant constraints on cosmological 
models (Bahcall & Fan 1998). 

In this paper, we present interferometric imaging of the Sunyaev-Zeldovich Effect (SZE) 
in these clusters, which provides a measure of the gas pressure integrated along the line of 
sight (Sunyaev & Zel'dovich 1972; Birkinshaw 1999). The change in the observed bright- 
ness temperature of the Cosmic Microwave Background (CMB) radiation that results from 
passage through the thermally ionized gas permeating a galaxy cluster is given by 



where Tqmb is the microwave background temperature (2.7 K); ot is the Thomson scattering 
cross section; and m e , n e , and T e are the electron mass, density, and temperature. The 
frequency dependence of the SZE is represented by j{y)\ in the Rayleigh- Jeans limit, f\u) = 
-2. We use the SZE to determine the mass of C1J1226.9+3332 and C1J0152. 7-1357 using 
the methods developed by Grego et al. 2001; in addition, we present interferometric SZE 
data on MS 1054.4—0321, a cluster of known temperature and mass at z — 0.83 (Gioia & 
Luppino 1994; Donahue et al. 1998; Hoekstra et al. 2000), which provides a standard against 
which the WARPS clusters can be compared. 

The SZE observations and data analysis are described in section II, and the conclusions 
drawn from these data will be found in section III. Throughout this Letter we parameterize 
the Hubble constant in terms of hi o, where H Q = 100h 100 km s _1 Mpc^ 1 , and we use the 
cosmological parameters ^m=0.3 and f2A=0.7 unless otherwise stated. Uncertainties are 
reported at the 68% confidence level. 
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2. Interferometric Imaging of the Sunyaev-Zel'dovich Effect 

2.1. Observations 

To image the Sunyaev-Zeldovich Effect in these distant clusters, we outfitted the Owens 
Valley Radio Observatory (OVRO) and Berkeley-Illinois-Maryland- Association (BIMA) mil- 
limeter interferometers with sensitive centimeter-wave receivers optimized for SZE measure- 
ments (Carlstrom et al. 1996). Our receivers use cryogenically cooled 26-36 GHz high electron 
mobility transistor (HEMT) amplifiers (Pospieszalski et al. 1995), with characteristic receiver 
temperatures of T rx ~ll-20 K at the 28.5 GHz frequency used for these observations. The 
cluster pointing centers and on-source integration times are given in Table 1. 

The interferometric measurements of C1J1226. 9+3332, and C1J0152.7— 1357 were made 
at the BIMA interferometer in 1998 and 2000 with nine 6.1 meter antennas in a closely 
packed configuration to maximize sensitivity to the SZE, with a 6.6' FWHM primary beam 
and baselines ranging from 0.6 to 14.3 kA ( 6-140 m). Typical system temperatures, scaled to 
above the atmosphere, are ~40-45 K in an 800 MHz band centered at 28.5 GHz. Observations 
of a bright phase calibrator were interleaved with cluster measurements every 25 minutes, 
and Mars was used for amplitude calibration (Rudy 1987; Grego et al. 2000). The MIRIAD 
software package (Sault et al. 1995) was used to calibrate and edit the visibility data and to 
output the reduced data in UVFITS format for subsequent analysis. 

The interferometric measurements of MS 1054.4—0321 were made at the OVRO mil- 
limeter array in June 1996 with six 10.4 meter antennas in a closely packed configuration, 
with a 4.2' FWHM primary beam and baselines ranging from 1.0 to 12.0 kA ( 10-120 m). 
Typical system temperatures, scaled to above the atmosphere, are ~45 K in two 1 GHz 
channels centered at 28.5 and 30.0 GHz (2 GHz total bandwidth). Observations of a bright 
phase calibrator were interleaved with cluster measurements every 24 minutes. The MMA 
software package (Scoville et al. 1993) was used to calibrate and edit the visibility data and 
to output the reduced data in UVFITS format for subsequent analysis. 

We flagged data from baselines when one of the telescopes was shadowed by another 
telescope in the array, cluster data that were not bracketed in time by phase calibrator data 
(mainly at the beginning or end of a track), data for which the phase calibrator indicated 
poor atmospheric coherence, and, rarely, data with spurious correlations. 
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2.2. Data Analysis 



In order to properly model the cluster, we must account for any point sources in the 
field. To identify these point sources, we used DIFMAP (Pearson et al. 1994) to produce a 
high resolution image, using only data from baselines longer than 20 meters. The resulting 
synthesized beam sizes, RMS noise levels, and point source detections are given in Table 2. 

We perform a quantitative analysis of the observed SZE profiles by fitting isothermal 
(3 models (Cavaliere & Fusco-Femiano 1976, 1981) and point source profiles to the inter- 
ferometric data directly in the Fourier (u,v) plane, where the noise characteristics and the 
spatial filtering of the interferometer are well understood. The spherical isothermal (3 model 
density is described by 



where the core radius r c and (3 are shape parameters, and n eo is the central electron number 
density. With this model, the SZ Effect temperature decrement is 



where 9 = t/Da, 9 c = t c /Da, Da is the angular diameter distance, and AT is the temper- 
ature decrement at zero projected radius. 

We determine the best-fit point source positions and fluxes, as well as the cluster cen- 
troid, using a simplex algorithm which minimizes the chi-squared statistic (Reese et al. 2000). 
We fix the cluster centroid and the point source positions and fluxes at their best-fit values, 
and calculate the chi-squared statistic over a large range of 9 C , (3, and ATo values. For a 
given electron temperature, T e , the (3 model then yields the gas density profile n e {9) at each 
(9 C) (3, AT ) point. The gas mass and the total cluster mass can be calculated directly from 
n e (9), assuming that the intracluster medium (ICM) is in hydrostatic equilibrium with the 
cluster potential. Following the methods outlined in Grego et al. 2001, we calculate the 
gas mass, total mass, and gas mass fraction over the (9 C ,(3, AT ) grid, and report the total 
mass values for which the x 2 statistic is within the 68% confidence interval (Ax 2 = 1)- We 
evaluate these quantities at an angular radius of 65", where our observational techniques 
best constrain the cluster gas mass fraction (Grego et al. 2001). For clusters at z ~ 0.8, 
this angular radius corresponds to a physical radius of 340 /i^oo kpc for an Qm—0.3, VIa=0.7 
cosmology. 

We can estimate the electron temperature directly from the SZE data by finding the 
range of T e values that yield a cluster gas mass fraction, f g , consistent with the mean value 
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measured by Grego et al. (2001) for a sample of eighteen clusters. To determine the gas mass 
fraction at r 500 for C1J1226.9+3332, C1J0152. 7-1357, and MS 1054.4-0321 , we scale the 
gas mass fractions measured at 65" to r 50 o using relations derived from numerical simulations 
(Evrard et al. 1996; Evrard 1997), as discussed in Grego et al. (2001). This calculation is 
repeated for a number of different temperatures ranging from 4-18 keV, and we report the 
T e 5z values that are consistent with a mean value of f g {r^o) = 0.081 (Grego et al. 2001) 
within the sample standard deviation of 0.04. To determine how sensitive the derived mass 
and temperature are to the adopted cosmological model, we repeated the calculations above 
for f2 M =0.3 and ^=0.0. The SZE derived mass and temperature are reduced by ~3%; the 
overall change is small because the effects of decreased distance are offset by a compensating 
change in the f g scaling relation. 



3. Results and Conclusions 

Synthesized images of the SZ Effect toward C1J1226.9+3332, C1J0152.7-1357, and 
MS 1054.4-0321 are shown in Fig. 1. The SZE decrement is detected with high signifi- 
cance in all of these distant clusters, and the locations of the SZE and X-ray centroids are 
consistent (Tables 1 and 3). Using the SZE data and the analysis techniques described in 
Section 2.2, we determine the temperature and mass of each cluster (Table 3). X-ray tem- 
perature measurements for C1J0152. 7-1357 (Delia Ceca et al. 2000) and MS 1054.4-0321 
(Donahue et al. 1998) are also shown in Table 3, and we find that these X-ray temperature 
measurements are consistent with the values inferred from the SZE data within the stated 
uncertainties. 

From the SZE data, we infer a total mass of > 2 x 10 14 /i7oo M Q within a radius of 65" 
(340 hiQ kpc) for each of the clusters shown in Table 3. These mass calculations can be 
checked against values derived from gravitational lensing, X-ray, and optical observations 
of MS 1054.4-0321 (Donahue et al. 1998; van Dokkum et al. 2000; Hoekstra et al. 2000). 
Hoekstra et al. (2000) infer a total mass of 5.4 ± 0.6 x 10 14 /i7 q M Q within an aperture of 
radius 94" . To compare the SZE and gravitational lensing results, we calculate the gas 
mass within a 94" radius aperture, for an electron temperature of 10.4 keV: M^f (< 94") = 
3.7±0.6 x 10 13 /i7oo M©> where the uncertainties reflect the statistical 68% confidence interval 
in modelling the SZE data. The total mass is estimated by scaling this value by the mean 
gas mass fraction (Section 2.2), from which we obtain an SZE estimate of the total mass 
within a 94"radius: M t s Q f al (< 94") = 4.6 ± 0.8 x I0 u h^ M . We find that the total mass 
estimated from the SZE data is consistent with the lensing measurements of Hoekstra et al. 
(2000). 
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Based on the SZE data we conclude that the newly-discovered clusters C1J1226. 9+3332 
and C1J0152. 7-1357 are highly massive, with a total mass of M? f al > 2 x 10 u h^ M 
within a radius of 65" (340 kpc). These values are comparable to the mass inferred 
from SZE imaging of the z = 0.83 cluster MS 1054.4— 0321, which has been confirmed by 
X-ray, optical, and gravitational lensing studies. These results demonstrate the ability to 
determine cluster temperatures and masses from SZE data without access to X-ray data, 
and illustrate the power of using deep SZE surveys to probe the distant universe. More 
precise measures of the temperature and mass of C1J1226. 9+3332 and C1J0152.7— 1357 will 
be possible with deep X-ray imaging and spectroscopy, which will be obtained within the 
coming year by the Chandra and XMM X-ray observatories; with these data in hand, the 
SZE measurements can be used to measure the distance to each cluster (Reese et al. 2000) 
and to further constrain their density, mass, and gas mass fraction (Grego et al. 2001). 
Additional, independent constraints on the mass distribution in C1J1226. 9+3332 will be 
obtained from a weak-lensing analysis of wide-field imaging data taken with HST and the 
Subaru 8.3m telescope. 
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Table 1. Cluster Observation Log 



Cluster Pointing Center (J2000) On-source Integ. 

Name redshift a 5 time (hours) 



C1J1226.9+3332 0.89 12 fc 26 m 58 s .O 33°32'45" a 41.6 
C1J0152.7-1357 0.83 01 A 52 m 43 s .O -13°57'29" b 27.8 
MS 1054.4-0321 0.83 10 A 56 m 59 s .5 -03°37'28" c 43.0 



a Pointing center coincident with X-ray cluster center from Ebeling et al. (2001). 

b Pointing center coincident with X-ray cluster center from Ebeling et al. (2000). 

c Pointing center coincident with optical cluster center from Gioia & Luppino (1994). The 
ROSAT/HRI X-ray center is at 10 h 56 ro 58 8 .6, -03°37'36" (Neumann & Arnaud 2000). 
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Table 2. Radio Point Sources 



Cluster 


Synth. Beam 


Aa 


A5 


Observed 


RMS noise 


Name 


(r uv > 2kX) 


(arcsec) 


(arcsec) 


Flux (mJy) a 


(mJy/beam) 


C1J1226.9+3332 


14.1"xl6.1" 


260.0 


-39.3 


1.73 


0.152 


CI J0152.7-1357 


13.7"x23.1" 








0.221 


MS 1054.4-0321 


17.4"x23.3" 


0.1 


1.4 


0.98 


0.092 






-161.0 


1.9 


0.63 


» 






-25.5 


-86.9 


0.35 





a Uncorrected for primary beam attenuation. 



Table 3. Cluster Properties derived from SZE measurements 



Cluster 
Name 


Cluster Position a 
Aa A5 


r-pXray 

(keV) 


rpSZ 

(keV) 


M to ^[<65"] 
(lO^-^o) 


C1J1226.9+3332 
CI J0152.7-1357 
MS 1054.4-0321 


0."2 
-l."8 
-7."8 


12."3 
-9."2 
-5."3 


6.51^ 


10.4l|° 


2.7 ±0.5 
2.1 ±0.7 
2.3 ±0.3 



a Offsets from radio pointing center (Table 1) 
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Fig. 1.— Synthesized images of the SZE decrement in C1J1226.9+3332, C1J0152. 7-1357, 
and MS 1054.4-0321. Left panel: Synthesized image of CI J1226.9+3332, obtained by apply- 
ing a Gaussian taper with a half-power radius of 1 kA to the u-v data, yielding a resolution of 
99.3" x 87.4" at position angle (p. a.) 32°. Contours are multiples of 290 /iJy beam" 1 (1.5<r), 
and the rms is 190 /iJy beam -1 . Center panel: Synthesized image of C1J0152.7— 1357, ob- 
tained by applying a Gaussian taper with a half-power radius of 1 kA to the u-v data, yielding 
a resolution of 151" x 87". 9 at p. a. 4°. Contours are multiples of 480 /iJy beam -1 (1.5a), and 
the rms is 320 /^Jy beam -1 . Right panel: Synthesized image of MS 1054.4— 0321, obtained 
by applying a Gaussian taper with a half-power radius of 2 kA to the u-v data, yielding a 
resolution of 73". 1 x 45". 5 at p. a. 2°. Contours are multiples of 120 /iJy beam -1 (1.5c), and 
the rms is 80 /iJy beam -1 . 




